{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<!--BOOK_INFORMATION-->\n",
    "<img align=\"left\" style=\"padding-right:10px;\" src=\"figures/PHydro-cover-small.png\">\n",
    "*This is the Jupyter notebook version of the [Python in Hydrology](http://www.greenteapress.com/pythonhydro/pythonhydro.html) by Sat Kumar Tomer.*\n",
    "*Source code is available at [code.google.com](https://code.google.com/archive/p/python-in-hydrology/source).*\n",
    "\n",
    "*The book is available under the [GNU Free Documentation License](http://www.gnu.org/copyleft/fdl.html). If you have comments, corrections or suggestions, please send email to satkumartomer@gmail.com.*"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<!--NAVIGATION-->\n",
    "< [Calling GDAL as External Command](06.09-Calling-GDAL-as-External-Command.ipynb) | [Contents](Index.ipynb) | [Unsupervised Classiﬁcation](06.10-Unsupervised-Classiﬁcation.ipynb) >"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 6.10 非监督分类"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.colorbar.Colorbar at 0x22bcf4a3160>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWYAAAD8CAYAAABErA6HAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzsvX+QZclVHvglqlYwQrQGWgPMlOQalhGEtfID4THgpcIoELUWtAJVa8BrYCmBFdFEFEsYBDaC2vWjHK5YbDkE2sCFqAVW6g2thJDpRsHIhgqFCaIijLySGEqDEGstpsWUZAY1iEKLQF0i94/7zr3fPfeczLzvvap+9ep9ER1d7/7IzPvr5Mkvv3MyxBixwAILLLDA7OCz7nUDFlhggQUWaGNhmBdYYIEFZgwLw7zAAgssMGNYGOYFFlhggRnDwjAvsMACC8wYFoZ5gQUWWGDGsDDMCyywwAJjIoTw2SGE/xRC+O0Qwu+EELaNY14TQvhgCOEwhPDuEMJKrtyFYV5ggQUWGB9/BeDrY4xfDuArALwshPA16pjfAvBojHEA4B0A/lWu0IVhXmCBBRYYE7HCJ0c/L43+RXXMf4gx/sXo528CeF6u3KWptnJMLN33OfFzwmV8+qFwr5uywGnhk8/AM49POpvP8pk/86PV9/Lpy0vAsz8DfPIZ1Y5nf8Y8bhx8+qGA5fs+Ye47+tT9yX0AsHzfJ6q/P/kM4Nmf6RyfLUOuCaivcfmBO9lzeb933NEfX8Ezj0/wwJf8mVsGHyv3lcs6+tT97edAWH7gTvsa+LnQdf3lHz318RjjA9lGJPDiv3c5/vmfdt9HC//vk5/6HQB/SZv2Yox78iOE8AwA7wPwCIB/E2N8T6K4VwP4d7k6Z8IwX3rO5+PvfsGrcLR2BQBwsnp8j1u0wGlg6eAylvfvtLbd3r505vUfrV3Byeoxlg4uA2jeN/mt2yiQtq4M77p13N6+hJ3BLXPf1uF6/Tcfs3W4jkeoLW8Z7ds6XDfL2jr8nmQdSweXsb15o1Pnw6pe+/zvwSN0PmO4u9G6f6my2m3/wk77+P5LXcPd1+At1G59j+T4J1//mtvJiyjAn//pCf7VrS8rOvaxR574yxjjo97+GONnAHxFCOF+ADdDCC+KMT6pjwsh/I8AHgXwdbk6Z8IwP/P4BEffccXcJw9xgfOPk9VjHOFKy/Dp56uN5WnVf3v1Uqserlf26Y6k5F1cOrgMDPz92pgNdzews3kjeY5Vhme0gcrQyX45tvQ72hncctsi98IzysPdjbruZNuxXpcDNEaYf6fa9djrs5dxTxBj/EQI4dcBvAxAyzCHEL4BwBaAr4sx/lWurJngmD99eQknq8c4WT3G8v6d+iMBFt7zvEM+9p3BLdOjnjZOVo9xtHYFK8O7nXcr93t5/07HW5ZR3tHalfpvy7B4xir1fotR9fbtXbuK4e5G0bmperYO1+t/Foa7G3U9R2tX3A7hZPUYw92N1vcr5eu2LR1c7pTDnUjOs58VhBAeGHnKCCHcB+AbAHxIHfNiAD8D4JtjjE+XlDsTHjPj9vYlLB3c61ZcPNyrkYnQAzJMLkEfr9q6LvGcgWP3mBJUba/Ksob+08Zwd6P2VrcO17G0Zrdbe8s7RI3wfoHer7G9eQPD3Y264xEjnbrmpYPL2ELTDi5fnh+3Q65N9kkdcn3tjuYJt957gAcBvHnEM38WgLfHGH8lhPDPAbw3xvhOAK8D8GwAvxhCAICPxBi/OVXobBhmNfmywNnjXo1MWkZ2HzV/mcLJ6jFWhndryiF3rLdde3bcJjnP4sB3BrdGnmH1exKjrHnUnKcoBk/arqmT2uihbZy5Lpu3Xk/SMGKcc9ea2i/3jamQ2lhT24FmJHWE5n0Q7vyxZAvOFjHGQwAvNrb/M/r7G/qWOxuGeYELhaO1hmeWySTZXtpB9J00FJpEn8f1l0AMSt05qPL0BBmfU4KUcZbORE9S7u1f7VyD5bF6GO5u1IbAokbE2OaejWfwebtVhp4MFMh7UtIBzxtmzjDzQ5qWF3eaE0oLlEHPwO/tX633Le/faXlGp90OizuuqA37t+Bk9bj28Ngo89B+e/MGhmh7lvUEXwbi3XrGWSbO5F726VBSkE7GK48pCDH2Hiw+WUN75jU9Q9dWqVtuYGt1HSvDuxiubdQc9oxRGaeCmTPMwPSN80wM0+cElpdYAj3jvoK7o/mE2bhH9dB57Uo9zyEGY+9a1YmIsdZcuAzztboAGHmMZJTlGDZ224ZETP7WRkwbMIFWSuSoCQst/prKr9uUMcqM1HFahdMaWSQUIVZnOa+YCVUGQ9QZ8vc8wOMyzyMm9dKWDi4nvbPTgqgxuB3awGoqRYyLnLe8f6dlwPkY3fGkVA7eeXpCbutwvfbSPQ20tMtSOFjgcph7Zv01t22cbzBnvK3JwNS9Ypy2amdWkPWYQwg/D+DlAJ6OMb5I7fshVDOOD8QYPx6qKcc3APgmAH8B4LtijO+ffrPPB+alYxGIR1P6EVkQRYSMiKxJHgs5z7pkdMWTSNYHroMy9q5dbfHhADpBUEJjrKhjLH55uLsBGIbfozpSBm5780btyUtbr998vFWXNREndESuLssDHxd8rufF86Qmt7t67g1nf3v1EvDKsZtyblBCZbwJwE8BaD3lEMLzAawB+Aht/kYALxj9+2oAPz36f4FzDvmgJvlABUxTNRM7aaMsxtGbCOrDU3tel3VtfCwbvhojIyNc6NHalTrIgjnZrcP1ynjvA1vbiocdQ9WxdbiOFTSaak0veeoIHVDjUQcex13aNu50ltDmkVPt3EKbY9/evAFstsucJVXGaSFLZcQYfwPAnxi7fgLAP0U7YccrANwYJfb4TQD3hxAenEpLF5gYk1Iq0zDK40BUDh79kQul1ri9falVlpQt5TCFkarXAncMTFPoey/qCm3sdNCIF0Ai1yEGuQ/3y6MPLn+aXjEHpej6c+07L8Elp4mxJv9CCN8M4CjG+NsjwbRgGcAf0u+nRts+ZpRxHcB1AFh64DnjNGMuMMsh59y2aU7GVp5l+fESyJFqmw71zqEJMoHJK+vtqYmnrcP1OiJw6eAyhgdtaoKDZ2qv2/BUh7sbWEb7Grz73TJeN5sAjc4+A62gDbQNqvz2ovL0Pul0lgAMDxrZXen7kpIHlui65xW9DXMI4VmoYr7/e2u3sc1M1TXKzrQHAPc98tD46bwmxL02jGdZd9+6tIZ0Wm21DG0p2NvTz66vWkRrmNkwrSgVgHiXHkUgFIbAMiqp9gnVweWwgiNnoDg/hi43dS6X75Uh4OAVvU23RY5L1T0OnXJRMI4q40sAfDGA3w4h/AGq3KLvDyF8ESoP+fl07PMAfHTSRp4mZtVbvRfQeSDOMvObBWv4r73iSeR2lkpjZXi3vg+ctEfAqosU2Khsb94ookO0cU9F6mmU0BCpY0rOF5020FAxfZQbuo4Ov2woRi4qehvmGOMHYoxfEGN8OMb4MCpj/JUxxv8K4J0ANkKFrwHwZzHGDo2xwL2BZehyx9xLeB+8JLoSozBumyVp1srwrkuDsDpBjh3ubrSMZUnypZJ25gycaKUto6XLZl010NAtNfWgJGryu0RmJ1ptuTdau8116wRJuu0sB1w6uIy9a1db5wpPnUroNI8okcu9FcBLADw3hPAUgGGM8eecw9+FSir3YVRyue+eUjsXOAVoQ9AnJFqjlBKaRkCJbicHH/B2UUmk6rq9fak1UuBczWKIt7YbDrmOSCMjIVI5TXV0tMeJYI9SuZ+b50L9FipGaIelg8u4vd2OsNN5N1LQHLM2lDrwpSRKUODlJpEyOSnTRUGJKuPbYowPxhgvxRifp43yyHP++OjvGGP83hjjl8QY/1aM8b2n1fAF+qPEqxsHZ5GuE2iCRHTkmIBTxqaS2VvQSg2GVZYYCTbYHLJtHZuCd+/FY8wFq7BHXhnhS9jevNEJZJHjhZaQv736tVfN99sagXFZpSMZ6XAsyshry7xjJkOy7zXOckJQ12X9tpLv9AWrBqbtFTeTaGXlTnJvzTwXI2NdKzP2+5UpUW9bGHnG+wCQ72h2BrfqHA5W6sycQbZUDhYkMEUHq2g+enhQUQj62JIMdqk2cApOkcJxeXoUxGoPbmeubsk+B1TPdWu7m7L0omDmDbOVE+C0cZYTglr5YP2eFGKUeZheeo0lkrmKDigL8JjmvZWosBROOzuZpdJIGVxNB8i2HEWh95sTZw7VwceUqjw4P7I4BpquqDukQVdvzfk1UkmQOBKQ84oMsYGdQfsaG9XME8m2zwNm3jDLh1ySX+FeS99mETwMH8dIaWPseYde+aeZpIi1yALOayGjjFI9LWBH9+lw5xxWhnfd42uvUCUmErCHqTPynRVaRn4AYLPZZ3Hl0oFY4dbmhPOgfZ5O2GRx5tIp7AxuLSL/ZgnC4aUwyQz9vYRItMSI6t+TlAu0I9j6TtIB4xlVaf+4iXC8MjV0+Serx7U8LTVhpIf6VhSeoI9RBvIyQ8/Qaq+zz72zOktRRwiEhkh9I30pAzle5w7R+/gfgzt60YrXeulM8qZ5xsx7zH2x8Ji76HtPcpF5PDQF2nRLKzFRD965BCXldZIGJY6xgjG8skoxrvHQBs1SPpTU3ZKlqdFNycfuJS7KBakAqGmOFpXjKFGEm+7otkfHi4GWd4n55ouAc+MxzzOu33y8/jdN8FC+dCRRIjGT4yQNpuQw1p7itDvJFN9uXR+n6fRgqQ4kNamnGS4N5pB8EbkUoIAfFl4aXKI7Sytsmu9FqcEvmXxjQ1xCu/CEoVeulCPfxUUyysDCMN9TiKyJIbKtcRLoWODycwElSweXi7XMuXZNO2owFRzDncQ4s/dbh+udobbQPlqFYLXFg7SplJYSGmYcSDu1gbeMshyfuld8zTnjrMsSaijVqZVI6i7y6HfuqIzzBO/F66geek7YaR65NYFqKCdkVZI+HwLL03jS72T1eOoqCDYs7NFbBkeWIUp5ypYygjE8aBs4bZxLJhItQyz5na0JVLMdJFObxGP07oWlkrAy4Xnt42M5zwhTFDr6sEU1HVTPaW/fn1wtWQB2HrEwzA7OUuFxGhpNSw8N3DEnX/oYUq3QOIuFMqWe0kAWndi+73O0dOVSt0wqMrRx3RncAm5W+8SwNMeMlqFSRkyXwUmNSiRunvFaOriM6zcbIykre/N+zQOzN5uK3mOlhFyr1RZWXVjl8ehLU0s5DfS8YkFlzCG82fxJaRFRWdxL6GsopXv0R62H2h4X3OqERtF0ubKt8/U2i/eXPBwCHsWkNNGl2znHRelz1G0CumqPEmgOXKBHanKd+riLZJSBhcfs4l4boGljGlLCyntrL79Uikn1zLmISC8fs9ZWcztaw+9EXYISZULJdg7UWN6/g63V9RG/3yT1KY0a1FSBtH940FV56Pbo1UQ4AKWiYuzRCXv4QFf9wlw3PyfZJhSIR9HoQJaLiIVhvkCYpLPJUTup1bP7RhoC3bb2oWA0fSGJiLQcy1uOqq9R9qCH5R3PcNBEuLEhKqnLCrOutw3a/G+u0zAj825y9F+aSpBRBBtoNqrWe5FK+VlLBdUyU81x8x/5t6AyFugFT/qVMsqlmCQZUu48K6TZm/DUx3ra3hT0JJrXMbGB4qRDHsZZronbJGk1S2BFe+o6eXsJTeK9DyVBL7ngmHlCiPGeLR5S475HHooPv+577nUzFshg3IRKPIFWoo8eB3o1a70vNbPPkjLJCSHbWwEayiiMoxawPM/SbbJd4Blkbqtck+70dHpNqw4JHy+9Tm6zrODNoxcdfMQGvztJardp6eAynnz9a94XY3y0qFEOHvlbz4r/6taXFR372CNPTFxfXyyoDAU9JF/k32gwrgqjpWrIJDri/MjyHHLGus6RvGrL1HIQWkG4Xr3unTbKtYHb7BSVhKwYrVUQHs/qQUfMMXRb63mBEW0CNO1PJVGqaY1EPVa7au4cd0yj3AQ9ddubokuYpnns9eYhc4UFlaGgvUHO8bvA+DhZPa7VE7n7qZd8ysE6RiIRPSVFqp06IIR5b6FAjtau9KI1eEHW0nMsWNfqKUk0ZLt+x73Q75xH7uUX4WctfLGoWnRbxGAvvrM2Fh5zBlaqzIUXPR5KlSEnq8fAPpLKAKBtwIe7G/VKIn1D2zsa4u1mmSOuQ/OpUq+mCHKctDWhxdkTPbWC5L/ge2gFcOg6ZRJP6449tQZfX6d+dCdopc5W4IzSK1v3pJtX5eIGlGgsDDPBMril3OhFwaT3QjxSna7TK9OjMyzeeGt1vCAEc2JrlO6SjbPJ744oAgk+WTqAmdJzC+u4vtnuMLi8ige222OV5ZUj8CLvOKDEM/5yr7e2HS960EzEcVs4EMbjiAXyHonMse44jE6Fr+eiYEFlFKBETnURkFq0tA+sgJCS2XpOlGS14zSCEIQKEcWEV4eMrLTG1wPTAUCTAMh7t1LURqpdDDbQqfLk+fCqN7pj5HuSa6tcK5eXa6OAtc8XCQuPmTCOwZ00cOI8QWb3J03l2Q0V745MRGUhS0UxPcFUA5CnLniCKed58bAfyE/MWcl+ShUWzZJcdtkeRC0B+FnpvG3ccTD1od9jGdGIxvoI3dEA16E120JlyPbtzRv15KrUo3OgSBi7BSt0fJ6x8JgnRJ/w1nnAaaTylIkf9sBkkg2wJ6tkfwmfLMeXDIenNWTOeXhynXrSKxVOzei7pJUuq85dMTLOetWUPqNEMexclvzmCUJ+DjWVMeqYvLB6OU46FC/H9rwha5hDCD8fQng6hPAkbXtdCOFDIYTDEMLNEML9tO9HQggfDiH8Xgjh759Ww+8lUkOxxcxyP1gfJHOVjZfexWnlT+gTsCGGRhsub+KLcf3m4zVVIkaHpWMcWGGhJCAlBzGgfC0CCQ/PgQ07n8vXoqE7JMvw64lVq555RYnH/CYAL1Pb9gG8KMY4APD/APgRAAghvBDAPwTw347O2Q0hPGNqrT0lWBzaadWzQAWttdXQtMCknpLFe+rfqeQ8Yrz4uGl0DFyGjB7EQC8ddHMs9wVfD1+DqCgsmVou9NoDlyMSOY8jvr19qbWQgR41aGrkoiHLMccYfyOE8LDa9mv08zcBfMvo71cAeFuM8a8A/JcQwocBfBWA/ziV1p4ytGRJoxVNZQRa5F6gi0R5aOh7I8EmnGfBkpGxlzRONByv9swKAuu3BvO4rRSUTipMy4iWGG9WQuhOSuduzgWeWHI4TykBAMO1jc799RQUJUhJ8Jhq2RncqvNmy8Rp7SFvt8uURE+nNUKaRUxj8u8fAfiF0d/LqAy14KnRtg5CCNcBXAeApQeeM4VmjA+W7chvj2tLJR23JnOmgVSCoPMCfR+r/MrN7+3NG53Vo/WEj2eUeTLMmiCSZ6mNby7EuqXThW/8U2UOD/K63J3BrVqaJvmapX4tdbM6hdrgFSpCGNw2z7Dzqta5iU2vjbJvef9Ofa0r+3fqOQJpd2vOgFbo5gnLRRKjDEIIWwBOALxFNhmHmck4Yox7McZHY4yPLl1+1iTNmAp4IoKheTBZ3847/yJ7xX2RmmAqlZyJAbAmYWW4Pi4m9dCkTalrESOcqsvK0yHbBeylChUguZRLJhelvYwSKRyDuXZ9HOfOOCvq8DxjbI85hPAqAC8H8NLYZEJ6CsDz6bDnAfjo+M07W3Dwg5Z0caat1PmngfPuLTN4JML8sf4tL6YVCcaRY7wMFHt6/JtzE6dgqRpcL9Joi4YXAWhFEMp2nTpT5yVubaNz965draMOLToiBx558P3SUYglig8uk0eSre9jv/pPRyXyubpc6Ww+1OvKzifGMswhhJcB+GEAXxdj/Ava9U4A/1cI4fUAHgLwAgD/aeJWniEWXu/poaJkuuHG2vDq39YiowBqA6U5UG2gUy85R5x5XCrz0bozSM0rsBH3qDF9L8SIMe+rs66x8a7rUBrgncEt7KHyUrksr521wYei6/aR1Bczh21B9MsrQ3texqJS5G+5bt5+tHYF+IDfnnlBiVzuragm774shPBUCOHVAH4KwOcC2A8hPBFCeCMAxBh/B8DbAXwQwL8H8L0xxs+cWuvHRGoYZX2YfaLdFsM0H9Pw/D0DUJqDI6U7rycJDQPGhrWvxlc6lpIOP5W+VMCUQcqD1fe7w38bC8xqSu/29qWaGsnRLfw/w5M9puZj+F5dREfpQuRj1t5KabQez5T3eTkWuTW60JNYXqIeSxZnHWvxrp43auVb0NFulpIhNzmYa0/JStG6DKBL9+TumT7fohNS12/tkwnnEp499czMXB2E3D21oiwX+ZhnFDlZm4aXmCiFar/BjRVgYZS7YKOTMjAcmt0yrAXhuPqZ6vX8SrXHltEeJ1cDXzMP0fW1aIPsGdCU4sEzysLd5q5by9xkKS4vmpBliCx5lL/3rl1tlvOi9vG1lWjTW9z86NxFPuYZRWoZo2mJ0T2VxkVE6T1NBSrk6B09rPYkYX3A70NpMhzmrCftkLl+/puvlYNJvDJku0xO62vILbmU6gSBspzjOlBF2qKXxrLa7v3W5Vq4iMElwDn1mAF/WXhGTv+b28/Ux0X0gkspH4YYYNajyghHEuKsDO9iuNb2lpb371QLkxpD2473jO6La2nPPaqCVR8CNpbaYx3ubrQUEPq8VvsKFSDaW8xxxnX5a9Xirdw+rQeukwNtm0V1Jjr1qNBTkGjv2UoLKhCuXAJDxHNPqVg07aEni5t7Nv865gvBMVsoNbYX1SiPAzEIet08oJm996gMLYtLUR+WUfPUAdpw6vOtOnNBFHrbNMKHPa5VB5uk5IPaiDJPzPc+x/PKudZq46W8s5TjSRglz3UqqlOPoqT+D73yx2aGYw4hPB/ADQBfBOCvAezFGN/gHPt3UAXg/Q8xxnek6jy3HvM0kDO6C6PsQ0+QSd4F/eHqvz2+0xoGV+UddwyOtw6dBev5efKxcQNKWucNJs/r4WHrsJKzpULTeXuthBiM5HPDq3WEK6/MY4FzMsux3uiyhLuvZX6D0XWstb8taacui0c/R2vApClnTwEnAH4wxvj+EMLnAnhfCGE/xvhBPmiUM+hfAvjVkkIvrGEumQBcGGUbrM3lD8fSKHvnWtCRfEe40hm29zWeOmiCO4BJIvtS51rvlqd7zr1jHeOFxohZkM5SOi8O2qgWuu1OaHNwiq5bgq1SShJXXWNMuu9du9oyyqw11+e0KCiePJ0hxBg/BuBjo7//PITwu6jSUHxQHfp9AP4tgL9TUu6FNcxA+qNYeMs2mB/mZE7T9kJbSYMcikHD8yL5JedgFYuK8OqYdC06i4MuKc+T7mmIZ1tpoCu+3oraqxYZ6LanoxxR/3s5P1IyRq0i4fr52nYGt2pvWcqwypZR1JOYHEefur/HhPITzw0hvJc27MUY9/RRo2RvLwbwHrV9GcA1AF+PQsN8LlUZp4mSHnnWeu1pQ6tb5G+WKeogCP6Y5J9WYuiZ/RR0JF+foB3L0HJQhnV9lpHWZVj1eMenOvVJl+fSBnJncAvXbz6O7c0buH7z8azEUwdvCEXhqZpKcn4AledtKTE0B69zLOf+FghVdg/wccnrM/pnGeVno/KIvz/GqG/+TwL44T7BdhfaY7Ywji563qCz5MnEi/wNANiXhDRobyccrV0ZK2eDhaO1K8Veq/bU2EMWD1EMDY+MUnplrVpYOriMJbVdTzRyR8CTeHoh2hS0t1zy7ukJzY7XrDoukb8JuAOWY7g85p31ZJ6mXuR3ndJzjd4bZ5kq6x7sDG7VipTHsnfgbBFCuITKKL8lxvhLxiGPAnhbCAEAngvgm0IIJzFGdwi4MMwLdFBzkeh6x23lRdlQXIOj67xjh7sbnaG1ZZQtY63lcRqtIf6g62GXUieSU1pydngKg85+dNvUSuajgz2w3qEV3Osh8H321Cp83fV5I7qBQ9S5PUJltdqy1m1Ha6J3gNGaf02bS7B37SpWcLfO3by1uo5ZksuFytr+HIDfjTGaoS8xxi+m498E4FdSRhm4AIZ5HC3uOMfPI8QAC4+8MrxLHl//xOkaOcMnHpeX1EgoA8sr9OooyftszS/wMTUnuroO4Lj+7UGMqy6vlahfnWN1Dp7M0DK+Xuen25N7BnrV763DdYCkezy6su6j/t2Hpx/ubmAZjRMwxAY4h/eM4GsBfCeAD4QQpMf4UQB/AwBijG8cp9CZ5pgn4XJT0VQLlMHj86z76VEAHjdsHS+crQytPU+QVQ/CZ/eBNjb6t74+vfacwKuXqYsU+Bq8dzQX1SeGTreFPVzv/FSHad0Dq7PQ0bFynBy7d+1qls4a7m7U//RzYbWItbbivUaM8SDGGGKMgxjjV4z+vSvG+EbLKMcYvyunYQZm1GOelkGdtYd4nlB9cMedbUsHdnIeL4+DcMPM42pNq6DUoxKPLxf40QfsdcrfYlzFKHvReiVGuOPFG5pdab/QRX3ywfD5uj5OnymdXf0MCJaUUDhloKG2xAvWlJGnIEk9I62a8ToY/j1rHPNpYCYN87gGVU+SnKbk7TzJ6cbp6IQ/5XMso1kbAoMisD5c3gd0M8C55SvuV65JdwalRtk6rjNJNhruX998fGRwqn3LoxBoMULcUcnEmE7WU5KZTu7Jyeoxrm82Idaaokhdp8c1Aw2NU1MoiYjJVj035Xc7Lap135nW2Tpcrzj4w6a+cTvOi4aZNMzjQjwbGYKfplG2DNcsguVJVqLyFJb37wD7MIMPrOGpJeOaNvTkVyrjGtCdVPPa5RkMa5t4spo7lomxjoLDuRaPP04N/bVDoL1P6zmtDO/WuTRqjh0+DcJ1ae5cv++p+vXEpXWu9q7LoiZnZ/LvtHDuDHMq8VCTjPv0jKV8GH0kT9OoT/7WnU8pLGqi9DxLfaE9rRIpm6Y+2MB6C9nKkFsvsQTkV5DmiSo2At7El14V2wv91tKyVjIfVZ+W1KXuSaftCnrlErkGTcEwdga3OsmGhNpIdVQ7gyqEWwJV9D7dLu9ZaOMs++XeLGDj3N2bnEE6bQ9WPtrTWhGbob2jk9Xj0erSl3vLjkrrYOghOUMbFL38UepYOZ4n1eReptQI+gPveGeqHouzZu9Z86KtQ+70AAAgAElEQVScV3h5/069QjO3WdpYG/o1YGfQ5nDrthXkkLauma9BOqWdwa2kPlmMY2mucu5QpXzrOeVGhpbixOSJ1XxCn+9HTwIv8jEvYOJk9WxWxPbKt/InWJCPedwOpG8I8tLB5U6+4JIAE76XoqbxhrTMw1rgGX6GZfSkHsuA63vmqUN0x1lC32jlgnXfGCXqpOX9O/Voso9qRu/X902/56V0Tw5eOy1cRE763HnMs4LSl+q0Jgl1mVqWxiOLOnJvv/uhsfSMt6dGJlaCGqENjlApNzRdoL1CPp4hnpTkZwbQmTQSCkR7utY16e2CncEtM7ERn2/RM5LL2VplxEvHCbSNK+cYYXpK8pBwJ5VqCxv2ozV1HYNKqjZc20jSKbozYXpJR/elkKJjUsZcj4K8Z3fRsDDMcwLOdyArUrQj56qPn42CwIuok6G4zgmsdb9A+0PzJGCyn8+RxPlioGWSSoyylRyJt4nB0XIwi5O26BQvAEPSd7LUr1W+4lE1bSEGvDZstF8CJSR4hnlcDqDgSEExzix3k7Z7xlMCM1LeskVdVIEzMDvfPgoR3sfvAD8/8xkN2tpxSbZ/kVCySvbPhxCeDiE8Sds+P4SwH0L4z6P/P2+0PYQQ/rcQwodDCIchhK88zcafB0zbW84l85FoPfGSV4Z3eyfNkQ9OtLzL+3eKlh8CyoadYlgEFtVS0mZNS7CB2hk0UYFyLRo56ZZOfKTL9+gWXhXbK5+vj41tjr6apjdp3ZPcHI7XIQu8e2K9t9LJW9cs2yrt/OWWrvwioIRjfhOAl6ltrwXw7hjjCwC8e/QbAL4RwAtG/64D+OnpNHOBleHd+qO0XnIxwpWKovrIS7ll/bGJMddGYHn/Tv1xWAoAD/LxWoZKG6O9a1exd+2q207LMHCbNJhH9spMbWODaXnce9eutqLW2Iu0Rg46itDqHIRn19fLAS2pa9aKDS4vd/1eu73yBdImnVdFILSIlzpUOjKuSxyMi0hvFC0tNcoz+isxxheNfv8egJfEGD8WQngQwK/HGL8shPAzo7/fqo9Lle8tLXWegjhOE+ytWUoG/q3VEVY4sTeh00R5dbXL3lJAJcNZoCzIgpUCvMpGSfmW0UnJ61LnW5JML3+GRW2k7klu0VmdXEjzrx5NlFJslLaN9/O8Qyl9Ic+tpqNoJW2v/TlYAUbTWFqqz3J206ivL8ZVZXyhGNvR/18w2r4M4A/puKdG2zoIIVwPIbw3hPDek+O/cCs6rz0l5+oYZz/D6pzEaIlCRM9y8wdwe/tSbWhSapJqEumKOZzNnaeHmaXetG6vnJeS6nlt0P94mLy9eSPp5Uv9SweXW3w9j1S4bOnAWnpdMmq5EYVnZDW83MV83cv7d1yvVv+WtllgI7i9eQPL+3fMEY53vryLjBSVVhrKzvf4onDN0578C8Y20yUfJZveA6reyzpGew59vGcdmNH3XF1/7lgBy7nY+KbKKm0fe7/sqemP3NSRAvWMP3umJ6vHWS8OaE8QsiqDJ9pa3taoXLmuFVoxuVWOuhZuR10+aY653Qxr1ez62mnizdPr8rXwRGY1UdrWrm9v3qgzy+nzrRVS9PVasCbCZNJWJ4hvqT1QdWIrwzvY2vYDTVL1ex1Ijm/Wi9pak6AC0UPzhLIob8bxpOcd4xrmPwohPEhUxtOj7U8BeD4d9zwAH52kgQzPCFrHWeqDUqTK1cekDGrJvpK2NpKqNn83PNhoHqDxQWwdrtcGEUAdlSb7lqCUESN+OgetZpAypF2sJuBOQSsEWMbVGm6LosIwqPo66mvNBDkILKmfBc9YpToMwE7wVHuiim7QHQwHuGxv3qjUCdho+H6dUhPUkW52uXHdDlE76DZbxxbpnnsYVFaIcCfnRQxedIxLZbwTwKtGf78KwC/T9o2ROuNrAPxZjl8uBQ+HlvfvZI2e1vFOCnmZcnVPGzwUFJqBaQygWY2DIS+58H2WLInv6e3tS73UG57HyZ4118nbmR7wyqulUsZEXO7jtThR/TdDPD8xGloKp9UffCyPFPTknz5HQ0+wWm2Vd9lTLuj25GDRExZK7nOuPk4fwGXp+YO+tNdFQNZjDiG8FcBLADw3hPAUgCGAHwfw9hDCqwF8BMC3jg5/F4BvAvBhAH8B4LsnbuDo4R+tkYSmR6J2q7xxDOu4OSpyOFnNJ54X3q42buwd0985j0Pvk6Wfao9le701Sac/divYwkuEIx55Z9uBLEPfzTBXGztUz1tWr7DuuTWZaD1bpi52Brc6q3+IURZPujZ2g25Y+BGawBe+tnqBWsP75+vsBpzcabVRB5ikvHu+TtFMW2HkfI0AcLKd1iaXQkZJ3vnD3Q1gze8ImL6S694ZjcQWXnOBYY4xfpuz66XGsRHA907aKIYV4SYf9lmimdQ4HW9ZKwF0h5R6WeWDl4Q6SwfdBDxeZ6TLFe9bDGRqAk4bB0E90WYsocTnyocp5y7pfavrrm7bGiV4SztZNAAHoIgHx1SKoOUBDyTYobknWq1Se996dWjVyfFcgRyzh7ZMMPXMpF2caznlNHjPkQORcvAoEOuY+lms5suWb0s6mFS2wIuCcxf5V7188yWhO1mtkhO1XsZEAhxruCwGTJKaa89JGyddloCNU2qiTNcvSYl0svWU8RQwVaMj+YAuv9ppa8YQt+pTE40nq1XuY+5Y9PVyeZomWcHd1m+T/7ZGHoZnyzIzgR72e9DzFDxKSHm2gFp8NbdElsFjM8TQc6eUg+aZF5jRJEanKZErpTFyE40e99gXfc8b7m7UwST8Egv/LBy0hgSHWMqHkjqZR/WQ8spsxUF1H6XdVnCCGKU+UV/WdbHEzTrOOkeG2d6+kvbJ8xoXwi/nDBzzuKnJvdQ1eHMVqU6LodOSeu9KKkHUwkBXmEmPeZLJtWkEpXhljKvASIE9HZ3XVkMHl0iyIB1goL0ell0d4QqGB349LU9LPnBKiLOyf6fFg1reodSplSs5z3e4tlGtzm0Elmxh3eSVLWjPTtQfR2hyS2vPt6OgGP3WBo/36dwcnmGUpEzyN9erQ735fqXu7Q5TJ46XnlNeCFL3NCWB07DmDMw2DdKG23svG3XSE2UNOseYScM8CaahmEjJ5cYpP+0VNzPzsmKI6FEZMpkmw93WPvgvOdBoXauJra6mtFOeUZZQSLLkEdDO5StgaVbd7oN28hquJzcxxMdYOZ/Fq8XNdhk8sQRUnZ9+dhYlxPVqeNF/lrxPX4tQFzz64BB775nrUZnIJlOJfaxr4M6EM//pyVDrnngdCB8vZVvrBqaQ6jAYLLHDB7LFnnvMJJUxKU6TCmF5UqlEqdSYizfFQ1/+wFLlCL0h/7Tmmr1Xi+ooGTrmPjSvDEuGp7fxsNoql6kPCx6VcLR2Zax81Jp+2Tpc79x/fV8FqfvEhkvO15N2ntxPjpXjtw7Xm0yCh+utd2VncKt+DzSVwtF/3G7pyLzc0BzhOS6N53V6uk5rodg+OZzPO+bOYxZMg9LQYN0o11EHiSRmxd227Lf3iSRr79pVc609ywCJvMzKWAZ0Dabl+VgTfa3OR7Xfol3Ya5R2rhhGebi7gZPt7rnaSLTaYowKWH5necMybB7ubnToDYu+4LbU1MOI+tHBIR36yIAndyv1JHMqiK3DdXNisPbI1y7XE8sCa66E/x6uieFte77iWVtRpuKBs6LHqks/YwmSalYzSUv55O/FKtnnEGwsT9M4y5BSlt0ZN8pQv5Csl/WgVQy1gVldzxo42WYNcz3elWkJgIaVuNKSWukoQwCda9k6XDfX0tMfrKdHrnNDbDeGl8vWxwPNS74zaEKmuRPxPD5WSXTKLpCY1QufFnK042h42VlILddkURR6H4BOdKiglVtage91kpNWz0qe8e3VS1gZtlUlqSXNLgLHPFdUhjXhBJR5KH2gPdNJ0Od8Gbp6YP2oHD9N8LBYJrQ0bcK/dwa3sp0Et509LctYSnkpakK3h98JLdHzaCJdx/bmjQ61Me13KkWLeTSNTNJKsiZRtlgaaa0i0vCiMU/LwdF1NQsp+7hI+ZjnxmO2Xh7u4fVQdpJ6JGJr3JWnS6BXpx7uVoqIZhKv7bXqj44nS1KKAW+bnpDiVJzA6KPfb2eBk+AWxs6gGzgh5QFdzrWu72DDpEtYAcGTUhwxKNdUtxPdTooNASsruA6G11nk4Ent+HlI291RAtoJn2SfPhbgRFVdvbjVDkCyCqIzkpD6uKP13h/+W99X3l4rdKy2H/iGl7ntBZUxB5AXaVKj3FJOjJBbQXhcWMM4MYSAHYXGQ0nZv7d/tUM5cACAlsbJ+VuH61hauwyg4QyFQ5f6qg/lTmspIh7G8hDeCwryDBGAbNRh514NSKGBhkIRnl72cdSoZTgAO7hG9ulIOc8LbSkvDPmi1CGdQL28ljqWOyIuWydAynHRukz5n71QntjW2zyYhlqt91jXdbDRosfquguCqcQpWVAZ5wwp4zipUeZZbZnl59n+SQIINLwPweoYrHbytfIMvgxJU+utWdvFy9F6Xh7uW9dvGVU+rk4TOQJ7RLmhc2pIzhQE/50zMFKnHmUxvyz7SzxnT1kBdK9dwCvEeGXpNrPCItUW/mdNBm9v3qifGXvI3tJPubaVQMr2vk/efhar0s8SZtZjPg1eaxxIuGqqLakkRNO4jq3D9Za6QSs2WGpmrVAt3gkPg+UcTRcwjy3aX8vIiCcni6hyPmCBXsj0aM2nUvhFtBQluUknPlePYqQdekVu2afr1OVpqV7L80tA7pGncRYsHVzG7Tq50LF7j3ib906xobPUNALJ8cHXJedrYy/3WXJZeKH91iRyDq0RiyqPHQm559OeN5lVzKxhngWjDKDF6Vptyhne1LlWORJE0dcTSalRtg7Xzegwi6/cGdwCbvoTdfJxyAdSfzjYgNU5WcoOHcFn8dLaIHu6Wp37gq+Lf7eMJLcto5iQIbiH3LOSekuO88r39lvb2JjtDa928mcATcIsnsPIwbp/XC7TMylYkkwpP3WfgZHscfPGheCY54bKOG1M0lGUGmWG5vmYlrCOYzVKrr4OH+3s5986X4L+QK3heW6NP6suhgRO8G/+V4pSoybIzf5biZhy4OCMUtQTg8a18r3xOqXb25dqeksWpK1GLu2wcK9uvV93knvXrnaor5J7wqM1fW1ayim/vUncecXMesxniXtJm+hhMg9D94btlJIAKx+aiTm98oX2YNijlEkwbeR5klQrMaoMau0w4BLj701CSRu5U/EMLatDPO1wPcSmdnnP1KJmWMHCZbKCwA1gKQQbltR5HmXRWYZr5L0OD9pBLHr0IMtO3d7u3gvP+Oqyme6SZ7U3rJQ2MqEq+3nSV1NZHYpqtG9leBd7uIplpFfyGVcZcx6xMMwopxs0WLKWi/yT4zXEIHTWihupKixoj63+cEYBDZI0nWf52+qC7gSiZ5QFWr7HH5jFK3vGp48xk/p1InqNVISZ7qDM9owULpJAqXONhiqiL9VkdQbWiMJVVwzQGbXobbyP5Wm5hRj0tQit0Hq3RwazrhfN+6Elj9yJNIqe9mri3IGIZw803wjTehcRs0FlfPIZ97oFAPrn2BC1gvyfKkc+El6aihf3tLyB6zcfNyVlug0CmbiTTsL68CXngAXmj3nFlj4z4kxp6GFpypjJpBXzoAxLjcJKAwF7ZuzhWe3RZenUo61JNNrW9z3R163P73NvS7bxyITnHjxY9AHQvH/VpHKbOpF6b29f6tBNYlSl87l+8/GkSkXoFYm01B07t31W5p5OG7NhmE8RpR9RapZ7Gm2QcsQo6uiyPh6mZWAk+QzQTS2pQ4qlTR6kHFF+jHMPWrpgZ+LOu2ZRDFjSRO9c5mK118kjEc1bl1yHeIF83nB3o/eq03Jen/vZ12OU+2514Cm+Wp/P90z+XhnerfOOAM17tnRwuUoLOxpl5UYTwrdze3iORJ/PEaaTyPPOE+aeyujzEUzbCMtvSeepPc/OhMbo4/ESxpRCl6t5yhUjdaROOsRyvL4yJUneU7IsVc1LjqADWawOyOJLgfYq0zUXq2gHGR5bkYdC8+ilpoA2nWHmwC4Aq2M0X26pL5j35zZ6kZzWnENNacGmfIAm3azQRXoUxvSFVQY/M15KimkNL8KTr91V3oyoN1lv8SIEmMyGYX72Z06lWH7xvPwK0xoa6Z6fy7+9ms4DwB+W1x4eVu84qgBOyK7bBTRGT3Iqy0QNGwz5+CbhVNMJaCpoYynt4Ug9kxe+aX/AfN+0obGwvH+nk9ual3ySUYl1H0pkdhotCmTQ9eplvyxiALRzPWs9r4A10mzwtw4bXbDZ0Ry2kwhVsA2mQL8XvK6fJZer+XtDIpm7V9xht9+JG3js9UXFnGtMRGWEEH4ghPA7IYQnQwhvDSF8dgjhi0MI7wkh/OcQwi+EEJ5ZWt5ZzbgK16uXjx8XnlHWdY4DMRA5iZhe9UM+ELnOVOck5wBd2du4+lsNXmLJO4c//NxzsThRgdBFFl8MqDXuEtBt6EuDcN2epIw9ZZE6shZde+bauLFXKuVweX1yUeeepWe0c3SX7NdUkMf5y33WOvuLQmMAQKgWth7jxBCWARwAeGGM8VMhhLcDeBeAbwLwSzHGt4UQ3gjgt2OMP50q675HHoqPvOKfnBqxrzkzwJeNjVM2cHaTElp2JMNV8Sq0Z2KpKxjWqiDjwlIbiNHgjiHlUbfkVCNYk2WWFM9SPqQUGwLL4OgJJ486sM7X7wR7xKnyvPuQei7cXu++6GtJvbMldeWOsdpizXFY36Ruh1XnY4888b4Y46NuIwpw3yMPxYdf9z1Fx37olT82cX19MSmVsQTgvhDCXQDPAvAxAF8P4NtH+98M4McAJA0zcHaGjR++NsrjSOZKj88Z8JLZc6Bt1ISCyOWC5rzCLLVbGd6th6EAeg/PNWQZJ6ZDTlaPaw+1ZCjLHyJ7gpb3yOAhPZ8jZdWgCESPpmBvTbxXixOW+8b0g9AJ7PVquJK4wv3msYOu8av3qzzIvAgAXyeX47VLe/u6A2H6wiuHZaYAsrmtrWc77xjbMMcYj0II/xrARwB8CsCvAXgfgE/EGE9Ghz0FYNk6P4RwHcB1AFh64DnjNqMIpcbzXkpxWELnIdeRMK9YAp0DN8cDpjKabR1KRro2dgbNqhh90fGoMh9wHy48R8+wAUopduQZrAzvmusxlsDzpkvOG2eU09BejVGW926ISgdfEl7NBlNzzKkOh0OvPQ+eRwKSwfAiURljG+YQwucBeAWALwbwCQC/COAbjUNNriTGuAdgDwCe86wHow4KAOZLs1h6LdYEngeZ5Gqf3/wthlrkZ8JHWsoQoEsZWB8+S+k0PKOdunZvUrN09WevPK+OFI3AHRN7zvoYrdDYGdzC1rbdabC33VHh0DG569FDe75f3nm63dKB3FbLem1hvVYN6XbnjGH9PmXe744yZnUd9WTjoAnPZhtQ0V9Ncq7hQfkE4nnHJFTGNwD4LzHGPwaAEMIvAfjvANwfQlgaec3PA/DRXEGffii4xqLV2Bkz2ONQHzmcrNpLS6WS8reNZNvQykewMmxyR1sGStDhdmnozxnuJFSakeIULbDxGB5sJL1UrWrQ+7DaXmJJ7onm4vUq5LzUlP7oreerw8mta9fQZWsPmVUp1rlAM1qQNln3h//WnU9dx9oGJAMhjww4ZzbfG0uaV3cWSh0kqpK+ayFuHXaXG5OQbz1/clEwiSrjIwC+JoTwrBBCAPBSAB8E8B8AfMvomFcB+OXJmtjAeiHPAimPpkTVYcq+EuDZdD1Z5/GIufJz0Xs8BN0ZtJMUWbI8iQazVuVOlV8KNgZSv3RWntGX90NCfLWXJomgSp5ZSrrI97ukLKucPtvYi7SO0RN8pW1xOWnnd4r7Fg9Xluzyyt86XO8of1LvBr//s+KQnQUm4ZjfE0J4B4D3AzgB8FuoqInHAbwthPAvRtt+bhoNnQTjetqi2shxjCnkPCn2jFteHgn6aw96RFsM1xqlhR5+6raKd5LSR0sb+RgOxNghb0sfVy2kebcjh0q1ScumpMxGq9p0Ei3vb7sa/uoyGLW3Z/DamuOUibqdwa3WUlZLB5c7i5J6nKmc17pmtPNEcM4S7XVa7bOUDWLQ+B7qZ+eqS2iEIdecuncppKgk6fS8dKKlUsO9a02OGKZ/FhxzIWKMQwBDtfn3AXzVJOVOA2wQxulpS6V0nODFQupllOATAbdZqAN7Fe7j2iDKxIh0EJrbtfTPeqiZkkEt79+pgy6kTI3UfcoNa0VG53GaJfJGPqbmUgsnjMRIrlDAiR6+p1QBWuXCEZOciIoNJec1toJIdLl8jRoehSFl95k463DfiePcTsSgmXh0wR249cy3DteB7ab9ml6ZRYQQfh7AywE8HWN8kXPMSwD8JIBLAD4eY/y6VJmzEfl3Cph02NNnodVp6aGtyaGc4Zf6gXYAidchsKG1+GT9cevVoku9HqD9MXmKj9KE7Ryam0IzwVX+PFJLg6UkehpWnZUX2eaGxdMWWAZNt8HbX88jHDbRd4KT1WNg30+o780z8MRbSn/saZZ1ObJA7MroXW5GitVveed4dAbgPEnk3gTgpwCYL0oI4X4AuwBeFmP8SAjhC3IFzq1hnhTjGnaPB0xNFOYmEdnjzSk25EPylsQ6whXTuIpiY2X/DvbQLF1lDXs5ci7Fw2qk5Hb6GsSTrT/8zXYeBvm4ObsZUx4pNQe3Ybi7geWRwRdKgpdPynG6YpQsGkG2i0dsjWa80QvDUtNosLFvtcVRi5R0+Kw15rZ4oxq9kjeHhIvqY2VYGeRmWStagzLZmmZUMmu5MmKMvxFCeDhxyLejCrr7yOj4p3Nlzn12uWlB0wGaIpDfEuZtRatpyKSZZWw5U5wslLm9eSMZYit15xaH9SYLOSWoziQmnKbodVOLdFpeGKsPBDnvW9QEljZW2mpdq8XhttQE9DffU16QNAVPGcF1WdBctDUfIAacJz13Bk2gjrxb8n7xXICXGtUCr26i7xW3S9LUyr3RVIk+T7LLubQYUL8/JVGn0j7pWCcdmY6J54YQ3kv/rvc8/0sBfF4I4ddDCO8LIWTlJQvD3AM6tFUgE4SiCJC/+Xj9si4dXG5NLOqXvE+Og9vbl3D95uO4fvPxmrMtSdrE3B/Xb738nFY0Z1C9SUCuN6UkyfGJubbKMfxR67aVcpa6Y0qdr6+bj5WOLJUqNJWASsDHME8t2ydVCck+ua/yf0oJo9sjx8r3IO+j/K3BnaTeXqo26Y1PPqPVyaX+oeKEH6V/ez1rWwLwtwFcBfD3AfwvIYQvTZ2wMMyF8KRMqR7cUm3Iw+byrKASzxsTD49fdv7IJLm59ug4n63+yETmBHQ/TG4j/58LkhBvm725PjPrOdWFlKUpBDYK0pmI9y/33io7JUvU1IRuC7fHO1cfD7Q7D8vzlxGQlCHPnJ+tGPpUZ+wh9e7uDG7VGelEPqlzUOtFEfR1ybnMrXs0X4mip8/8xozhKQD/Psb4/8UYPw7gNwB8eeqEBcfcA5NMKGqOU/PG3kciOWhTGca0wdMGQz56/ZLXHKjBQ8qElBy3vXkjmY+XIw6tcqxVtC1DXeLRsrROpG362uW6eEVouRdHuNIK2ADaqUe9Zbm4foEse7WEJmeGvv+pwAjmxUVlM8RGHcwjbdnevNFZSsozcoxUZ2gpKzR0rmyvA5Pz67Ya9Xcm9xJgVdIc4JcB/FQIYQnAMwF8NYCfSJ0wdna5aaJPpqd5QImu2tI45yRTOW6zkwAJXe41BT25w+ByvOvTden1AXWCez6P2+0pAfR16PI7k4RiqI0ES6m1C3kizkohylI3d2JM3SNPJaEVJrmJ4tQzTBlp6x3h89jg6+yBlgadryHVLus9lGu0ZJ5bh+tTyfZ23xc9Pz7yHa8pOvbJ178mWV8I4a0AXgLguQD+CJWE+BIAxBjfODrmnwD4bgB/DeBnY4w/mapzQWXcI5SoMJg+0MO94e4G9q5dxd61qzVNIfu9mX0Lusxcmy1tdB9OWENPQHFZ1nUxr+rVxctAeUa5BCvDu/WSSQLhguUeWPSBDPm1CkNDlkyy8lTL31y+N6lbck21nM7g3b32WeV3FC2ZpaQsqsP6LcfyNXJdsxz1F2P8thjjgzHGSzHG58UYfy7G+EYxyqNjXhdjfGGM8UU5owzMGZWR8ybuJTqSt31/gq/FP9PQkF/m7c0brVW0vcgvoHrILOa3hq7WEJi3t8qnpDOAMYQftBPj8wrJrG0FGgmV3IuWB7a93pJTCS2yd+2qSf3otmqpn8XxHq3xkkX57G5Wuxo6p/08tXa7Y2goAVXpxJdQG6nRiB5V8L1jb5/fkcrzb9fFbbBWQeHj9N+6k9GUzsroHeDAHkDUIg19s4URFTdakWZl/w4+1GnB/GGuDPMkRrmvUbcMmzUMB5oAEACt1YAFNa85mgRkXi0XkqqH4aZnOCCPSRkrabcXPOKuwjyaXBOwJnbrsEpKw0vc120aANhUQ2cKeNESRCuSbGntMjgRD9AYwdTwWU/QLaG9UrkFSwve8WxH1yTQea9rLl/dX5b7ddpmUBtMZ/G98p4p8+1e+5ty2pPPomkHqvmPKjCkfW8seZ7UXZc76J4jy5sNdzc634J0lq0y0DyjrdV14JWdaucOCypjAuSG1QJrpWfGyaqftAjwZ6O10SiZtbYmsrwP2AIbBJ2nw6ojxY2nAm6sgArr+mRYLsd6owYul/nekmfi1e0h18Fb90MoHW+/ZVxFYwyg1q/3zcDmecBWLhZGKhOeQIyrXFdOeSTPUdQePBKxnu88Y6485kkwrrdcIumxAgmAtjeqX0ALWpKkPXSXyji43HrQO4NulJh4bDkDxGoNmYWvs4oZORlSRrlEecFtE3jqB65Te2xe5ymcuRzfysfsrFriTVjqtsrkIqteNIfuwXKhtPoAACAASURBVBrBaA+05W2OlC/yHolaRLdV3yugMgIyWhNFiI7sFKrDpm7aASlWRCF75dK2TupPZ0SYe1fmEQvDPCZKlBXeObxCc20Qt9c7Q0ULzOPK+d6Qlg2tluRZuYd1Mh09SZSie6zcImKsO203wDxoDpZR1kbL8t7YqGrZmRhi7qD0M2aeVura3rxR5y0G7DzVrNxoy/TS4fWdjnL1uNMRtO4BKPfIZterteYqWttXq3eQlyKr67nZ1HOyelyvtF7XTXI5nRMcaBJ2Nc+gfd+sTnZH5QC5SMZ5QWWMidJVRgTiZehJInlxdwbplYb5WK9czovML3Fq4kaXxQZbhpNWG4Q35TBuhhgFztc8bqLzvh9kyuuXsiQSz6NytFwxB1ZX8D3jFbtZhcFUgeSJ9igMPXHsTdTqa7D28/VZ4HfQ4t1ztJTWeHMwjLzjJZw3P5NZndA/Tcylx3zW6gz5iEvj+E9W26G5MlmUohL0QqdSb+OliFdofzwl98M6pjVMN5LxAFV9e8OrOMKVVtY0Kym9XPeOMWGpvXOgaxwkvzQPfbW3yzSEVpfodssxLaqEJmO9CV2e4NPPQKeqFA9aG3tv0lFTJqzpFeUEc/B6FFXfJ1WmZ1QtKVvqWC2Bs44d7m60lrDy3uuSd1VTMbOWxOg0MJeG+bSNsgzLgOOkUe5I5Gg7R5pZtAKDvSu+NlFw8PI7uQVLBfzxZxUlo2GsV/7W4TpWYCcTEmPBSXb0vTI/7IONTi7o1GKx0nbPWHhLHjEPLEZd+ForUEQUJ3w+X6cnFwPQys8sSzlZ12TdZ26LppssuG0f/Z1SmniTrKwukeM1hVSfOyVeWHcSy/t38OTEpc4+FlTGhCjpBMRw8rDWg3hIJ6vHrRwWemZ6Z1Dls9BllSoruO36GjiAxDKi+sOVtjC0sZF7UJrPQRuevWtXi6iFnEKGwTyynoiznmupMiP1DPg9KC3XUtKMC2tUkmuL5UHrd4ihA5Asj9w7t+4cDXi02TxiLj3ms4b3ssiLd3u14Q8lsKJePmet6wnqY6WsetKKcjjIMlOpCRJpnzn8HWG4u1HREGSsvQ/EUoKw0kMbZVmle2fQBAzkciYD1STU3rWr2MPVVvKkPp5YTWWoa9XbUue2ylH32dLqyvGas7coEx3sUXJtPDnJuaO5XqB5/3K5nPW1WvWVHCfHynugaavWMmWFeTMs2uyx7FnnHwvDfIrw+LOdwa1arSDJdRg8K623L0F5PaP/tUyJDUltHCjQxBra18ZYKQ+8a9OqDUm8wzw4J+zngAHPoOlO6vrNx5M5NXSb5Bp5mzXpqCffvElObdBqvlepMjzDl5IhluQj5nKkXXXHMOgaSqFBWsY0E1SkeWB5P5jbzsk4uT65J8vqfZDlvoC24kKX701gLlQZC0wF1seoZ/7lfx6y6+OkLB526+Eil8f1WryqUCWsa93evNE6/2T1uLUopgVRDLAiAejmgBDILL03wvBm4PU1bG/eyC4GkML25g1cv/l4a3JPOg4+Rkd2piD3gZ+pDolOjUBKtkm7Sugaz7P13hmGvB+sHrHeN68+eX+EtpGydH19lS8XCRN5zKO1rH4WwIsARAD/CMDvAfgFAA8D+AMA/yDG+KcTtfKUMI4WuW/5e8PGuNUTdIfrtWB/79rVammjfWAPtiHUnkVbhWGoHhwaQu8bHmyYw2teFTs14SZtkJwPFn9c8+SjY7huC6Ve0dHaFexdu9oJ6dXgla+1p629NBkBaPlcJ3iFPFBWcHCQRklQhB7Oe1yuNurciWiagSmWDuWC9gK91jPg8oRCYhWQPralehn4lIy820yvdML1E7hI3jIwOZXxBlQJoL8lhPBMAM8C8KMA3h1j/PEQwmsBvBbAD09Yz6ngLNQbnKhGsHRwGbe3K1mVqCr0fr3kPOdGFtkah0Zbs/ACzbOyofc+UIFWRujy9RJYul7r+gTycabCtlv3gMpdyUygdtQLKhkOgDrCTQyYJwMT3tS6r1JGrbJBN6LSgpzLRtzznHXH7N0rDhKywIE0qck5QbMuX9OZWkE34iHzby1TlHUVGy+5eX4lK3n34eDnAWMb5hDCZQB/D8B3AUCM8dMAPh1CeAWq3KQA8GYAv46cYf7kM8ZtxszDnIGnDoEjosSAaWMo51kvpSW1YrQ84JEB7obUNhBvTIb3snjrzuBWRW1sN20Z7m7gZFu8qZK70UZqGJuTxomEz4OeeBNuWzzAlKGrjRZ1jnxfxAACFQfOOuwl5J+JGNpK53vD5fPZ28154PJesXSvFC1DnjjXk+jljKW840yNmG3IzGk0dT+RrG8eMInH/N8A+GMA/0cI4csBvA/APwbwhTHGjwFAjPFj3lLdowUNrwPA0gPPmaAZZ4vTCl5hikK/vGxcU1IjDd0JAE3OAzbKllKDvfOV4d16WCvpLjsTcgftZYes9umOQK+UPA1YKgzLs2VNcef81SZ/hp4T4M5NJ4xfQjrvdc2zrrUNrZ6E0+ctwY6OK1k4lr1tC5oOaxldgyoDbMmkBT0RySojoP2el3jNUuZClZE/9ysBfF+M8T0hhDegoi2KMFrQcA+oVjCpCz3jqL1Zgly3tWS8nvnXK1to6Nl29sg5uAXwPZ4m9WOFZTWrrtvIE111GQcbnYQ4fLzbdpq1LzFAgvqaiQdeOqhSSepyRI3B+RrqNg1IQoiK9ri93eZaWXZXB8+s2rpv8dh5mSur3Qwx9gw2tNw5lhpLhqY2zBwcTiInnYAIIM+beOQd9a61OjNKC8udvlVfn3dgHjCJYX4KwFMxxveMfr8DlWH+oxDCgyNv+UEAT/du1DkyzjkDOQ60QkBLo8TQ9PE0rSHkieMZ8v7rm4+3vL0Uv2nmToZthC0jrjHOsLx1LsHN+1BQvqx3uHTQKBZak2SkXrHuA1NHuWfGeaU5X3RrEpDK37t2FbiZDnn2eGfdVqvt+jnpyU/ZBqAj38vl3WCkvnemkS4KxpbLxRj/K4A/DCF82WjTSwF8EMA7AbxqtO1VqBYiLAYnPZlF6LbJsk+ngZxEKQfJPVzxmd2IM64ntW1ncCt5PqATm/vlMoSn1qsvT/sD1EEc+m8LW4frrXzHqXeSKRpZDVz+yXWJbEyQ44t5YjdXt14BXOvLU+3Xz0wvo2W1NduZ9kAqGlTrz3nJsHnHpKqM7wPwlpEi4/dRLTb4WQDeHkJ4NYCPAPjWCeuYaZxmJ8LSIqAdTAI0CgAdtl0PbdEdBmqeT8AeCa8ULYEoTIskdbfK4yrxlFg/PW2jLIZnnLSRsuxVHaU4aOgDK8hHzgHa3DvQVWlYUs0WnTDo8viMet+mzfkzLH5a6k6pNFjlYfHpcm1evdoD1+8YX6vuNKXOi2KINSYyzDHGJwBYq8e+dJJy++A80R6TgoX6gK9qYE9HS+GYJ5aPUSfIARr5m3gpJdRDDvr8EhnYJHXKfZAPPDXBZknltrDeyS+sE0fV2+pOsWtMrLzUbAw1ry2diNx/i26xIhJL7gdTER4vbY162FinEnDx9eh3j+tbGd1Hzk9dH09RidwWvezaPGMR+TcHkLBiLxBA0JmkgR8Ky0N4Cc6wyhoXvJirt1/aNEldWlVhbfeOSYFpHa331cbEao/IEV1PkyBUT2fewSgXSPPN8r++ZqEQNZVo/RaMMyEnckMpi/NTc9IuaZe1/SLg3OfKmHVvOefRl3j8tadsBKvopXfqiKpRIIXWRLeCIYzAEo4MY6+u41UbbS7xbmujoToR8ZjqzoXOmQa9IaOH5AogolwZ3TuZfGOlhaTr5PstqKMMYed0ttqk2+BJ1AB7so4TR/Ud1bD00lJj5HhwrzNoqXIMNYZMpu4MboxW7hkFoVhrAqp3/mT1GHh9ydWdb5x7wzzrmFbHIR4aR9pVZR+3JF9a3qZRGrzC/KPn1enZeStrXG4Cidfvy2luS5FSmVhtAZoOUoJ9akOioh6tMvXQX7bVyzwVwqIAvE4EgKt44eP6TNTxs7L2pcqyOnCvk2E5JKtMrFXVt7bXW/TdaU20zxoWVMY9Rl/DLTkZ9PCSFRyioNAz3juDW63cElyGljmlhsYyvF06uGwmKuLjc8mGdF3TUOWUeIzcRqZNxvHOrTaf5kjutGRjORVQilLS7xLTFd69sK4jla+7NJf3PGDhMeP8TSBabeXgFPbSUiuaeMN5zR3WqxoDnSRIrEKQY6R+GYrKkJvDlZOpO9FNsNR3EtCb0OPfqTUI9cSWaLnFONQRgGjz8VxPX/RZE5E9TYtWqJU5PSgHnk/wwsrlHVkCOoE5W1ivg5JkCSwPnUAkypvC7dSh9AAWVMasQzjBSQM8xjXKlkEvmYSZBJa8SqMl5TIoB0+ZwOsK6pzIHlrcM9pc4s7gFnDT1sbqMGPen5oYGwe1sgDKKA18Y2hNbPHq0XxMSqpWiq3DJjGTBJTUbe8RZCPGWEve9GSflUoVQDIJUt0e9Xd75e/uvEeJZhyr7RW2WalSKVrOj+M0DZx7KmPWhjdnqbv0OGQZkuqkMVaOZKBps1AfpYEQ+jhvlp6Hs6nhspZmMfoavZTaoQ91YilZrOvsy31a18jUk9wnCcdPnaeh1SESeCLQOmyrXPmtKZ/tzRv1P6CZdJTOVOexlvZY5elAmGkvHXaeca49Zj3xci/qtzBN46y9Yv4tNALvO1qjdqkcyHoixfL0UrPxncm5jEFjeZaU7R0jHrcVeNFXR52a/PPUEVbbBJZXzcEPVj4Lgc5BYU2Gtu6DoYRJecz1MzFCqXkSkUdOXqcldFPn96js0pGlJ3GUfOB8nFxDQxM1UZq6nTJCvgiLsZ5rwwzMvlxuUqSuzxLctzurroZ2uNZwt8tK5A/4Q2cOTmh5YpoeoHKsqDHPsGpjxAmEhM8cN4IPSEenSdn1ckrodoi5QApWu2iFCd9TTwnRKjOT/lIbfgss19Pt51FAH+dGvxuW523BkuVprnu41l2eixN31VLGC4LZoDLmKB/zaXvw/PJqbr0krwPQHtpqRYeGjpgTyoKHs8BogpCWqprGaIYDC2SoXAJt/DikXMNaqkmuy/LeUxN044yU+ihBtEHm9mkPP7cSO5ejoakKhpXTJPf+6H0WbeIl2WKIIukiYGY85pQy4ryoJqbdTl7I1Cvfq+9ktbs6CntRlc72uDWMt4b0OaNRe40jj6e1fJWEGGcoCPnYU9ei2zKNUG0534qY5A7L8uRabRvd66WDxnteQTWRCrQ70FT0Xu7aNN0g9dfPgNrlGTBWbHBAj5TvRfml2mcdZ9EQ0oFoNcYSgKM1YOkArcl8Hnmw1n2hyjgrPPsz58Lw5nAa18B0hRiLEo9Yjj/ClfqD1fmYBTuDWzV1sLJ/p5a95YIXuJ4dWhmZI+UkanG4tpHkpD1Yy05NYpSBhIHnHMyk1Wa4S3hRIh5NL+n8Gam2W9cm26z8EaVgOVzd+SbmF1Levydx9NrP9NTy/p3W6jM6UdERqsRRejTII7aLkCh/NqiMDObBaI8LfkFzPNvO4JZ5jM41YA0lZUb9+s3H66EsTxQKpSH/OL+B7LfaJkEtqaAVT83hyfRSvOqk2N680VK2WO1OGUdtlDTNkDOqfTuclJLFO1Y85pSBzQWbWOfkfp+sHnc8eU+d4nHxZ6l6upeYDY95hPNCWZwVeJaaU3tqTSfft+3NG60lk8RrljIqGkPlyXA8Y9HTLh1cHq3g0U7cw+VVaBtmS5mht+vrsNQFFrwy+nrT1oQiXxuXK/fBypHM9WpqhJ+PrKK9rMpJtTl3PUIfWd/OuDy/lDVJZ2cF+YwLWd5sBXfxoYlKOh+YKcO8MMptMHUhORxK6IxqyaR2OUA1TNTLSgEJpQBzxdgAdwhi8LkdzGlrGRfzupacilfw5t8aLPHTRnh780a9okcprLK4M+OhP18zGxzmeFse9qCrLJAFWCVowjO6VvRdH0+9PscIGLEke155uYCTTn2s+1b8NdDN1eJJ67gsiSC9SJgpw7xAF3oyxpI4eVRA6oOqc/0q6MkWPj7Fq0r7hF/WQ2yRzfHkk9S1PArhbV2HIxmzjBW36fb2pZYxzPG5go7Ol+rn8ix9tZxjhb+z4ed7akXf6fP0tcmEnXWeZ2ytd6DPyNQ636vLnPiVqECM3qHV7uoquk1WtOpFw7ngmBfo4mT12FxNW0dS8T9L0THpIpdsaKz2rQzvtj469pBkQmvckHpPEysTjqWoV4hxOjLPEOr6SxUsYqRLj2V451nKDL2Ar5yf4vsZVsctTgHLI722WQoPqdeLQgXsCMCLJJUDFh7zPcM4fHo9sTca1glH6X1APPyzJlQ0dFCE9lq84JNWIAVoFn6UD9oKZNHepiQtkt+lsIbCO4NbrfSjfH+sYBcdNON565aXK3Iv4ak1LaIDWFKaauvaUrA86i2sd0Y9qVWwpY1aOgd0jQOXp6kFjuqr3yNVHy+R1iQ56n4D1jyKtO326iXglfb9mCcsDPM9wjh8OnPNqVWvG5RHSrFR9pYO4glDDpNe2e9GmEmINXs5nJ2tEzWWkObpduaUAHobt8k9f9DcA2/JJGubFc2mjbs8r2mspq6VE5raqDuPQdcTtjod3akzxaTBnbUsGOBRX9JRcRs4SZM+3lq04SJjQWXMGVKRgRpihHiSZZwOw1IppDoP69gSpCa/Smf9S9UPk6gRLAOeCvhIIZfASEZM8i+XSIqP9a5R6CCLmhhH5SHUB2Ab5QW6mNhjDiE8A8B7ARzFGF8eQvhiAG8D8PkA3g/gO2OMn560ngXY2FUGr6E1Rp7ISHXRChQYLd1jraQhH5onO5MlkgQi9QLQmnGXD03oCP54vTXv2IgMDza6k2/UDsv71bBUEt4xHrSCgtsg3ilfrxzHq7CkYHmsHnXC6pTWKukqV4WmKvT16vuhU3TqtuhnoPOdeBOcQFuVoj1vHilZk4apiEMdAXsRMA0q4x8D+F0A8qT/JYCfiDG+LYTwRgCvBvDTU6jn3KMvr6yPzyVx0asxA+mkOGJMNXQejfqjVHlx+cPbQrNckID1wHXUmSGHqzXaIymallzJ75Xh3awULqV7Tu3njmLrcL1DD0gbLa7eojPqc416eCKMI986fPGoHUwXSB1ajibttOYFzERTRi5q6zqs90dP5tYYtOki6TS4M7aMcivBksGHb2037924uuzzhomojBDC8wBcBfCzo98BwNcDeMfokDcDF1z3MkXoAA9reNxHjWBBGxY9Wcb/+CPReSbkH6/ezWAjI6t8m4txEvpytFYYNG8TI5KDKEdK6vfuiT6m1Ump37qd0gY+37pX0vnxvxz4WTH0qEXLIuUdKEFKp3yRvOA+CDHG8U8O4R0A/lcAnwvghwB8F4DfjDE+Mtr/fAD/Lsb4IuPc6wCuA8DSA8/524/8zA+M3Y4UTjua8KyjFfVEUipMm707PVxkwb+EPrMHLHWU6oD1cNPzIHV5e9eu1vWY+aEN45qCToikz7HCoy3u1KMmvICSjhJCebm649K/9THeNehnlErrKWW2FkiFT/tYHHWJJrw0HNwKzPHmDbw0pQDw5Otf874Y46PZShO474ueHx/5jtcUHTuN+vpibCojhPByAE/HGN8XQniJbDYONS1/jHEPwB4A3PfIQ+P3DvcQdXAE2kP20wQHceTAq1brgBKmJ+TjWxnexe3tRpZUYpRrI7XWJMnp0BGqTsb1m4/XgQgrowRI/NH2wdZhRaek2mrxrTLkrq7fvm6+ng7vSoEyOoGT5mpTEjamlbQxteC9b/oamVrg9up7pcvy7meqPn2+V7bcRw5GEWijzB3cyerxzGWXCyG8DMAbADwDwM/GGH9c7f8bqNiD+0fHvDbG+K5UmZNQGV8L4JtDCH+AarLv6wH8JID7Qwhi8J8H4KMT1DExTttQisdyLycnhNawqA1r4ky3UwJB+kB7caLPZW+Ql7Zi6sJDFUrejhjsM2T2jHnpsP5o7YpZl+ZGrWOEtknVw0ZY7pWFkoUFShU0qWtPPQu9T34LrdG30wTyUYnWslRAeZ7xe4GR+OHfAPhGAC8E8G0hhBeqw/5nAG+PMb4YwD8EsJsrd2yPOcb4IwB+ZNS4lwD4oRjjd4QQfhHAt6Ay1q8C8Mvj1jHrGPdFGZf+KMnPzPkqLEOdknJpA72Hq/V+Tw2hlwuSco5wpQ4iqFUbuxumB+Z5lh7Yo7Y0yR2jQaoBbfQkzSevds3Yu3YVWKOyR162ZzyHBxut69SBNLzNmkTU11iXq8LCtfdrKUb4HK/DkvzRMlLwqKOlg8vYG17F0lpb6+5BUxXe88zl4j4H+CoAH44x/j4AhBDeBuAVAD5Ix0Q04ojnoMBZPY0Akx8G8LYQwr8A8FsAfu4U6ji38F5S5js9sLoBaAygd46oNFIeJw/JJetZSTm14TACGST4gOvYwrqZj5kj5wQpAwXYMjB9TSZHPOgaZ87fYQWglI6EtNrEoj5k6O5J7FKSP52rRCsfeJRSkk+blRas5GHjzKOW7c0bGMKXBZpRiMQfa3WIztLH90knvpoBPDeE8F76vTeiYgFgGcAf0r6nAHy1Ov/HAPxaCOH7AHwOgG/IVTiVAJMY46/HGF8++vv3Y4xfFWN8JMb4rTHGv+pT1rzLYTyjLKHLHqWgJ/lS0jkdbZd6wTsReFSGXrU5VY4XrMBG0jMSkv/YgzYUgD+5yLSKVc7y/h0z2s2rzzOeHnhZKv0s5R7q/NilSNV9strOezJpEIfmpMdpmxh0oV24THmveJumsUrrGQfPPD7p5JLx/gH4eIzxUfq3R0WVzKt9G4A3xRifB+CbAPyfIYSk7Z25kOxzPKTJwnvRT1aPW/kvWLMpHrJ8yBWdwZ5cs589atEQyzmi/9WTduwpcgg1T1JZKTAtyqDkA7a8QYtyKQljtuq3FAxspKQuHkJrD3zrcL2dKa0Qcqy8w3Lfdga3asqg2nZnxGf7dEBqRFCXaSgmUpF1XObetaudY+sOi75BpnvGhTXisXJzCDw99ozahqcAPJ9+W/NqrwbwMgCIMf7HEMJnA3gugKe9QmfOME8bs/RAUx9bs85ZE3CgjbIgRV20qA1KNCMfYv3AaSjLqoDOKtWr7Wxk+hpqKsSgNHRbtdSM92v978nqcb3MkCcJWxnedVc5abUNfji3ZQR4n8XdejxshzumeyJqGsnFbBllvj+6bqnLUrtYho/L0wEpW4frwDZabdg6rKimWlNNuS4svpw7eOtvqz26Lg+WVE6onOq9ecI99x7g/wbwglHE8xGqyb1vV8d8BMBLAbwphPA3AXw2gD9OFTr3hnnWkPIqmWeTj9Ea8urORowYZ3IDUEvfBGLkSlJ9iryOPWqgG52mDbyVatK6bo72a5UHtLjgrdV1LB10uWCgktpZXrL2XEvQ4UeJt2WvGmjf/1Rnm1M9eAZV168liPq+l5YvsJwVvvcCa2JXyrZkkTpxkdWZW0acDbWmzKx3fZbkcjHGkxDC/wTgV1FJ4X4+xvg7IYR/DuC9McZ3AvhBAP97COEHUNEc3xUzASRzn8RoVrzlEuRka1o9IbmOge56ataHlzPKXvJ8mbAqjSSzroH5Q/nwSvLralnZcHej4y1q+ZbViViQzkEb9hwsvS+XM9zdwN61q8n7xfs8LzMVjWdBnp13DWL0uFPJzT+0qC4ymrpt/NubrJVnp2k0DXmvRW6Zm9+414gxvivG+KUxxi+JMe6Mtv2zkVFGjPGDMcavjTF+eYzxK2KMv5Yrc+ExnzFK1Be8T+eb0HI4y3P2PuiSTsozOhgZXD2Dr9Gk9kzXZakdUhndRL1hBV6Ip24FJgj4RZd6qxSqKrLMmZwdHthSPwvyHCwaRQdV1CobVkkQJSJ0hjcq4G2yLiN7p3JNUk/Oi2WwAbWWIrOMpWdAeRRgJTti6HddB+5cBCwM8ykg9eLqcFzG7e1LHW9TuMmlg4YzrYb3/bnzkoms1oeiNdK08KpV93B3A8sFOaA5YpKhaRvpBOR45pO1MQMMzfCgaZeeEGO5nxidjjZZGctSlCqLUiMUj4qojSRdm6hQeHEAfQ6AzjyA5el7I6oS+qvFZ4/q6+w3Mgjyb2+NQYnKvCiYeyrjrDGpnMcyzsD0KBkZlloSNIFQF/I3UJYcqeTjtfJBcBSbyPSkThnO6o/SqqtP9CI/p3HurWeoU15/iYRNc7BWhJ8ncbQoGevaLInapNDX5bXROlbz+xxFqqWfFwULj9lACd0wCYSH9TyAo7VKjdBFlQ95ZZST2Rve5bwLPdT1Fvjk4zloQeqwPDQdXCKQkQJLybShGq6RKmStCXbRw3O+Bi5fe+DAyItXQ2FtMIQy0EmJBEsHl6t7TpLDVq4HA55B9O6LNYmWk8vJ+Vr6p8FlW3RMKvpOrtVSSnjta63ETp49gNZ9zKEe0TiBT/OMhWE2MMkyQKVDXqt87hDEkGgvgV9Qi87gc2upHRkU7Z2Ivlaiumo4nVI97F3bMCVf3oSOtF0MsgzjW7mCOZ+z4tP5mgGYw2SL2sBIIy4dCdDw4MLdb2HdDMmW+ytRb6WqCG3E+HirU2iVNUDnvtTHweaXOcucPl7O8YI3vKx0fF4qqZCGxX3L8Z7DkDL0W4fVs2nykJyfyfxJMDNUxixF/E26Gm/fa5EX0ooI4wRFnKioJFRYz7577eLoPq0ntsoDymgLhnyUerLSU3tYqoStw/X6/mhKhlGyLJSUI7P+Wm8t7eO2liB3X1jhkvK25Rihc1LPsK8TUaqwsZ6BlgoyPIWFePL6u+JrslbdZqNutWWesfCYDaRexnHOL8HO4FadNChXL9MA9VDd+Tj19tyS8x5SQSZcLpfHwTFWpJueDNIeWkOhjP9B8ohD2sz3maPftCEQpLTS8lt7lq1rVMeLJjulNdZBQs2aed1jxevWHVWtg6Z6AbQmAdlb9SYdWZu+sn8H1zcfb52r/5ZztBzRM+geDSX7cvK6CEZ9fwAAG8RJREFUecTMGGatVpil3vG02yIftiga5D5obtTbXuoxWR2OBy+AYOngcrUOnQqp3aGAFLkmndvDXPZK/1YUAVMeVXu7XKM2JlKXHO+NgJgP15SDVYfAWuJKG+WOQkQhxQVr5EZwLK9rYdBESOp3RHcs+jkLdIczXNtoGWqvPVhtVnTnsrzzmGvXMj1rxe15x8xQGUDX27rXyE3AlICjtyzIUBWwl4tij6+O7JuQamFozlmgQ6U5cESUAhzgoiHGQLeVh9ElH9okH+PJapUjWef5HWfoL5AheapdOhijBB5NYT1rrSnOdbBenml9nL5OUUgwWCVihYX3uWb97llG2/Oy5x0z4zELZslTLoEeruvORasRNGQSyqMsKrQDTuTYeqZ7ChpPy7uT8nnCj4eW8puPrSmN7SYJEvbRLF81WikE++08z+IpWUPpmmoYJVRKpTvdOlwf5YFwdLKjD5w12VKHB01tmIEsg/bv2gN0lBsdQ2zI4vi9SRlfzY+3rjWhm5c2aCPgGXGdT3sL6x3ts27z8v4dYNNtelG9jFlw2M4CE635Ny3c98hD8eHXfc+9boaJ3BBKvEjPMMrLmTKcOQ47R/GU1JGq0zMcol5oyZ+oHdwh8JA3FVWmDZe+JkuLa0Fzwt65zBmzMWwNz40AFOt+WG3U16uDVTi6z+Kg9TXqZ211lLKPz2Fj6N1PRiqJlDeBuXftaidplPd8gbZahOvxpHf6vfHa/tgjT0y8Bt9znvVg/LuPvLro2F/9wM6Zr/k3U1TGLGLSIVSJ/rJUXZHa34fe0EPQ1Oy6N5TVXnqJMeBj9YSQbpdukxdc4HmUMtPPuZ4nGS5LeZ6CxMP25o1iBYSUI9SL7mi2N2+08j3LP534X/Z7kMAVrlt+e53H1uF6cukyvnfyjliOgr7/POnI6hNdv65j3jFzVMZ5Q/OypqmKPrA84EmNN6OV2W00BN+7VikVWitrr9nhws2+sjqt8GlpM3/Y3mgAEDWCnQNET0ZtHa6PgnC6EO225QFrlUGL30cjr9vbv9q6T3I8jwR4iJ97Nha371ExQDuUm4+3DJflhdb37aDNVcukn+e58qSg9pQt2saSHaY6MR3yz9cz3K0Cqyot+hNuGfOChWGeEEXeLPJJfRiltETfCMX6A9ter73JITZGCfKNnMOJcnPLWuljhxjx1YZ6YengspuUPiXTk98iP6vbjXaH2FGEGN6oFxUn4Mx+XpL7mspRyXo8mgewA0FSEYVcl2RrE4PJyhmuz/M0PfpG1BA5ekfKl7L6QiswvLzYQNnIc56woDLOAH1eqtJ8AOKdpbxMDf5AOSeFRLd5BkNg0Rol7bQS3Usd1sw/70u1yeJ6BRzOff3m4y3PzqI9rPqv33wc128+XvPsQhlZ180ea0oiJgbbAlMuJdC0hdAm3FEx5WG12foto5dUO/jZWM9K4L37ekSSAr9DFwULw9wD48wI516qUokU85upycLcRCJ7OMItjiMJLPlQ5AMXw6/bMAm0ceMoOZm09Np4snrcmUgr6WT4fulhtuzPXZvUZR1X2g59jtStf1v3aLi70dquPWpPtuch1TmKtNCC1CHtSU32Af3nUc47xqYyQgjPB3ADwBcB+GtUK8e+IYTw+QB+AcDDAP4AwD+IMf7p5E299xhnuNb3HM8jq7LO+fSBpjVKohVZogaMnxnPkuvlPrRUe1L7vQnCnVE6VGlLZ9kkKOUI2gbELNOYoGTUVAW1XYbk2msWjj3nhYqh9FQc00IrFzSagKGW7t7KR2KgE21I57lKndV27vAWpWFEBtbSu33gQ2Ne83nCJB7zCYAfjDH+TQBfA+B7QwgvBPBaAO+OMb4AwLtHvxdwoD/6lDHVuSaAtkHmffyRaXjDTkbOc5Jhswz5dXm5Yar20nYGt+qgDM/rZM24UC8l1yLtLPWKrWALLou5baZIpgFPuuhBB+t4wTsyauFgG0uaJ8dYuaJzowHvHnu8Or+jOeWFzLuk1nicJ4xtmGOMH4sxvn/0958D+F0AywBeAeDNo8PeDCQY/QVa6PtxS4IYD+yNiMHRBtPyloe7G/XkoKg1xkHuQ7aMMLdPG+8Ux9sXJZI+D/o5lbSnJLGSV2ZOiyxt4iCflB7aqsubUyh5hqnO0dK067r5WA8l0YvzhKmoMkIIDwN4MYD3APjCGOPHgMp4hxC+YBp1TIKSYf29RmWQ2lRFKrAktYSTp9Y4Wa1SYOZm3Lc3b2BvP22Qc4ZGlAt6WK6DQsQL5vzLK8O7GK5ttHJy1G0d5A2+NQwvCV6ovU60pXga3qhF0wOTIGeE+hh6XZZuX45eAdq0F79fmgay8oiYHHbme+Q2VkqY5OFzh4kNcwjh2QD+LYDvjzEehxBKz7sO4DoALD3wnEmbkcSsG2WAIwe7qSb10Dql8qgN/NoVs0MSGVluDTVRIngGILVPQ+RyrOvVWlotEeSVsdnQ5nhoSUzEx+kJOynPM9as77U6ME6sw/ejxqD/JBrD030Ltg7Xx/5wrWemIxctg1vX63DnQJMNMEVFlUxwsspk63C9WJY5T5hIlRFCuITKKL8lxvhLo81/FEJ4cLT/QQBPW+fGGPdijI/GGB9duvysSZpxYWDla2Y0wRhXWhF6QhkIh1iikZYPgT+IlCFOebEWN67r8WbchbfUeTv61M9g71b4Yl1uLhKwRCOsMU3DIvfEMnJWPaJ8SMn0eD9z/Z2JPQV9HI8Y+J5q2i1HrVi/L1Lk39i5MkLlGr8ZwJ/EGL+ftr8OwJ0Y44+HEF4L4PNjjP80VdakuTJyVMV5oDI0OOOcHvYLLL6Vc0MAPm3B3o81cy4fKSdokuAKzgNhKQf2rl1t5Uhgr1S3j+uyriWVf8O6Xv3x8rDbMixeus4S6oPr0/dbv3N9+NG+qgytJJnEWxdYdVu5vD3Kx3om1n3V0PfRwode+WNznytjEirjawF8J4APhBCeGG37UQA/DuDtIYRXA/gIgG+drIl5lETfnTfIoqxMOwhHLLAoCYnIEznYinGMfLi1kR9l/+KoM01VaKPMZQEqrHkb2Bl08zNLnRwKbqE+Z5S9zAvoAOAu2STX00TqdZPyW8aSDRtH1Hnt1WXI8HsSqZvmrBsJoG3cNe1gcfqTgLl3D3qiUN8/0W7nrqEu4xx+s9PEJKqMgxhjiDEOYoxfMfr3rhjjnRjjS2OMLxj9/yfTbPBFQ2rBVm/f8v6dVp5n68P1yu7zEecS9LBxkuG30Bolw1IeEqe8LDlWc/FcXx9Y1zTJZFsfWEEhst7duEP5cQJXvHJ021JIRViWSABF5niR1BiCRa6MHhAlhEcfTBtadVFFEaLO32yJ8VNJkySfLnuOuNlNJ9nxqKlcnsST40VFAaBOHsQrXgu8QA4pR66HaRStW01F1+ngGjbmpdSElCNGAUAnYKSkjGmhvi+Z/Blax1zzxaAUpAf9JHuC1nNDW1VjjWh4/87gVr3ALLczBWt+Q7AyvHshAkzm3jBPy2hWL7xt8KZtlOXj0lRFPTxdbQ+1dVv06tI13UAqBx2M0sjv6GMaLRsvxlEbZYFE2vXxKpke6PDII1XDEa6YHDmveF1qHLUqI2Xodga3OkZMX7e+/xZ3Pim4s8qVm+tEgH4fuza4KQWKrktPjlrUmIYXhGI5HvhAjws5p5i7XBmT8mne+dMMbkhBeGXW9Ar4JdV0hYdm2fcG+hq8j55XpNboGzChYQUv6HbwEJzrYy+a1QB8nuzzJgMnpRt05KH3XkxCP4wTVWgdqyd5+5Ql/0qvI6WZzn07Wu88ThnzgrlcweQ8qjAEnqHllUI8ZYZlgHkyK7VgqCgpdFm5ulLyPWuykJHyyrQH5SlA5Fhdhi67NBJOUwL6eEteV1KGIPUMWJ3Sp74+9U5jMlDQoivQfWa6rj4UoPf+bW/euBArmMwllWHJlc4DUnSJBTGUtdFNKB0kJzF/OEwncFlVsMfIqK+ut7j1jhpi1U/faHHSer8ED1hD9txvvg72ks0hMQ3xmYbQ5aUoAdmvOWupt1POoO3pe+lNU9F0Vhs9FNEZUzLKrbIoQVNO98z7rJzeqferyev9BOYd55bKyL1g580oA13vgFNX8od7e/sSbm9fqo8vGUZbARyecL/vMlXSHgte24a7Gx1jr4fMOujD6gC2Dpvcwey56WtLKQRSKDGUbFi1iiA1YmBvWBIH9QmX7gOtWBkXesVxwB6xeMfL/dELuAJlucgvSsL8mTDMy/d9ovc559HwloADM+QfG2fmNSXs2vp45WORhERiYIa7G9i7drUVSMFr4zGvKwZHdw5S7t61q/WHIu3R8AyLNQyWNugkPPre8DWWGBsv30XK6OWojj7lpYbrVnl6myhV+vC8p/F96Gcm2+SZadpFH6+flUgZpYyjtSv1aEQiEbWxvig5meeSyjivEA/yCFeq3Mtq5Q+WkgHp1J5CP+jEQNrj2Bncai01ZeUlsFbJ5nK8sr2PSKsvlg4u18tbecdrasqa9LNgcZxeZ6ZRKwmIIhIaRKsypuHZbh2um0tslVAwcr5st5Qlk4KvHaDOiPKO6M6DpY8po8rveOt3giqbZ8yEYT761P14+F434gxhJZYXw8OTe7UqgtarE2ORM0gscdPleVnpADtPgRWNxRN7leGs6hzubrhRgoyWd5XgxgUpTbSVVIjRWnx2DMjE49bher0uoBi9EoPsGUntZcuIRt8PnjMp1WJb0JPBJWD1jJzb6qiJX+b/pY3yPl3ffLzextCTnHqSWuZd5H3aOlwHXt/nqs8nZoLKeOZH41QnJVI4q3q8ulkvzNstQyvHaVlcylPOISX50+3qTBSikaBp3lB+b2/eSHpGmo+dFrhM7SULVy8yxL7D/JPV41Z4uWzrQ3ek5IEWdZGDNwk6rjxPYMkOvevU2mSrbs25WwZc4E1YcvTmWQT1zAJmwmM+S1iKjWkGoaTKEZpCS86scyR9peTG0LkShHbQHo1AJyEC0PY6UOZ1iZdU57ig9vJHyStMb2/eqHjs1W75nodsXUNtxBWFYbVZjKcFae/t1Uu9Pc7ae1NtXjq4jOFBXsKnjZKejJT3saZFlMFn2kdGBvp8/q3phj5gasIKDa/LHbTP0YE7faCNOauH9LvSdOoLVcaZ4NMPhYm8QEHp+ZoLndZESR++U2amLepAJvXk+FRODGtYqvlXUU3IR8Rly2+dU4K9LzEUy/t3WoEVvAyRVb94+VuH6/WEoxX0wXVouZWOYLS8Qk0BWc/B89w0SjxOLxBCl221Vasj2KBa7bKG+jLcFxUHG2l5FtZzyYHbpRUi/Nvy/D2lhQX2mvm56O9yUu//PGMmDLNgUgNZen5JPuJxwR+afjlZwZA6l4/NtVUm7TwwrSDDfD6evXztket2lqy3JioP/i3XoSmZ3NBepIEa3pDdS1jUd/ibi2rUFI9nnC26w+qI5V+qs/AMXe6dl/2l34aeXO4b4Zm6hlyHxkFU3Amzbnxe1VgaF47KOG20gzQao+JpcNPn2npRKxClm0jozogGqf7XBo4VGJqCsSYnrVWwLWOxYqg+6nNG6g9NuVhRaXoyT9fHw2Y9MSblSgck9wXIhyXLRF1u+S0e4UlbSgJDdgbNKtqTpAa1lrHy6AZuS8nEH4ecW9REpyyjLlauiFc83N2oJ0FX9u9ga7utIgHkHWneXZ2ydmdwC48lWz8fWBjmDCwjlcPJ6nHrhbIMqYecgJ49Xh2t550rnu7W4XrLcHqRV5IgiKHXchMDppPma3mTlnDp9eD4mJRCgz/2jsE0zmk/g+Pez3F5/06tJNBtra+FuFjLQHnwDJ7mvXNrCFrbvWM9FYw21FZnUTsE6rd3PNdjccgARnQdOvehotUazT0w3cRQ5wUzRWXMKixvt08kVZ/hV0rRUKrn5KhBLs8636JXNB8rsDx862/OFT0NpYBXZwq6/alE8xolXjUfW9Km3D2wPFmrXKY++iD1rpaWxxON404uAl3Kid+3i0JV5HAhPeaceqJUpdHX4IpXO4StQ5Z6JcBEUyFCb0haz+V9oiuoHl12rdAwvHZZzVpgqRX2rl2tDZunZmDFQR1luNZcJwcgWAEUpcEanqrCKtcKeLBgKShub1/C3rWrZj5os02DdFImoL3fomn0c6uv1Rkp9A0g0fdZAns4X7PH+9dt7aG6sJ6JaOhl1OBNyFoyvEnpn/OEuTfM1hA2t0p0e/jYNmYpiRtDT2QtKwoht0q11MVehgf2Ti3OWNcvyAWBSPtX0Pa0SwJcNHJGlD2wVLCIZUTlt2Ws+V7wB96hb5zJqD6BHB1db2IlbStKT9Mhmn6wRizSdk+Zw8aMQ905FLoPdNRf6v6kOsZ6cs941hzIY2u8F3K5cw+RnvWFZ3jGTaLClMJZxPt713x7+1JNcXhee0l5EizD8BQI3qSYdXxph5cDGys9QSa5rFNt0WWVQkvGUudaXG2flT68+qWs3HuvVQ4evcAopT207FH+aa/f61A8XJTw7Lk3zKmIuj7nitSLJ9xS8PS6/3975xpq11HF8d+fXFtR+7APJdSah1YhiNggWrDtF19J1V6rRVJEA1aKxIBFBCMXTCr4IYoKYrFEGmxKtcFH8X6pVlpR+qHRNOZpmuamRkx7TfCBFXwluvywZ+6dM3f2Po/c2WfbOz843H3m7jOzZu29155Ze63Zvt7U773BCFOoV2w927ctL1v8e5g/8b0x9hfjnLshmjLGccPhhRCPuEI91iVUxDqJy7wR8fI0jZRDI9rvAg77GRufUBepdvzvfXhhE6njnIoHj+vx9Ye/9bJ6PYf68Pv5xYKadBn23/8mDntLuQ18f/u5ksIoi7pjEbu24pDAWL64/rrvS4UXvGFOTesHnYb3uyhHTYiJR20x8TKa8UOXuvjeEL9Sm08KiUPSVmw9u8AQp+SEKqrDP+QK6/JZhHUXvb8w42y3YYkNW5OfO9wnxRffOJ/Y05Qc4vGJFU31Nd0ovDGKzzmvN38cQuPo2021A+lFmOrWeq7DG3d/PsUxxikXQl15qJ9wv/Dm0s+4pvQXl921aVfWHIQukc0wS1on6ZikGUlbRqlj1Ke/o7aVIo5s8IzyZDzEG19vnCcev7jxzdfAgpA4v2+djPHLTO/atGt+FTm3f5yA4pde9H1MJZXE09SUcU/7BucJp+xxVlmd66OJuP7U8UyN5gd1GfSbQofGOZ6VxG2HDGtomm4qKZrO0dSMBuZnYOEN11+LKT3EbYQDiKb2hhn4tGUHRqGfrZN0oaTd7v97JK3sV2cWwyxpGXA3sB5YA9wmac2w9SxG6Ix3G4zimwqnYItxpw59u34U8ew7L59bLznu7yAuF19f7Fqo+3243nOdD9GP4lIXkp+++v81jbj9Pt4d46fpO255D8CChYH6PWhLTeHDaXqdsaq7qOuMc5Nrool45NxUf+hvTfl6mzLuQldC3drV4U3CG/NUxEjsVkudR748Tpqqa7Ppf352MHWwdyGsuA9hglCXGdDW3Q78xcxeC3wN2N6v3lxRGW8BZszsGQBJDwKTwG+GrSh1Nx7WYNeFoMX0S0KoRrjpNYsHkcEz8fjFc6FO4XrJ8QVx7vpqEaN4VFqnA9/P0EedWkjItze3hvKmhfL6N1H3hKfVvCHbj7DCxJRUm/7iC9+qnXrPXczUwfczkajTl6d+G0anNEWe9PQvcBXE4Xth9mDY/vkQth1GT/gojH5tNV28dbMOH3o3lxjDwrWQfXLOiq29CwqF54PfL775eNn9DSNcPzruZyiPryO+4c9lYPqIl+7FOQ9i6yaBbW77+8A3JMkaXria5WWskm4F1pnZx933jwBvNbPNwT53AHe4r28ADi+6IMNzBfDHIgNQ5OiaDFDk8KwwsyvPpwJJP6bqxyC8GPhn8H2Hme1w9Qxi6w67fU657yfcPrU6zDViVqKs5w7gOuY7t7ftt9Cm6IIcXZChyNE9GYoci4uZrVukqvraugH36SHXw79TwNXB91cBz2Vqq1AoFMbFILZubh9JE8AlwJ+bKs1lmH8FXCNplaQLgA3AdKa2CoVCYVwMYuumgY1u+1bgsSb/MmRyZZjZOUmbgZ8Ay4CdZnak4Sc7csgxAl2QowsyQJEjpAsyQJGjc9TZOklfAPaa2TRwL3C/pBmqkfKGfvVmefhXKBQKhdF5wWf+FQqFwv8bxTAXCoVCxxi7YV6M1O0R2rxa0s8kHZV0RNKnXPk2Sc9K2u8+N7Ugy0lJh1x7e13ZZZJ+Kum4+/vyzDK8PujzfknPS7oztz4k7ZR0xsV5+rJk31XxdXeeHJS0NrMcX5b0lGvrIUmXuvKVkv4R6OSezHLUHgNJn3P6OCbp3Rll2B20f1LSfleeTRdLHjMb24fKWX4CWA1cABwA1rTQ7nJgrdu+CHiaKp1yG/CZlnVwErgiKvsSsMVtbwG2t3xM/gCsyK0P4EZgLXC4X9+Bm4CHqWJCrwP2ZJbjXcCE294eyLEy3K8FfSSPgTtfDwAXAqvcdbQshwzR/78CfD63Lpb6Z9wj5rl0RjP7N+DTGbNiZrNmts9t/w04ClyVu90hmATuc9v3wRCvjTh/3g6cMLPf5W7IzH7BwnjOur5PArus4gngUknLc8lhZo+Y2Tn39Qmq+NSs1OijjkngQTP7l5n9Fpihup6yySBJwIeA755vO4Vmxm2YrwJ+H3w/RcsG0q30dC2wxxVtdtPXnbldCA4DHpH0pEtTB3ilmc1CdRMBXtGCHJ4N9F54beujru/jPFc+RjVa96yS9GtJP5d0Qwvtp47BOPRxA3DazI4HZW3rYkkwbsM8dKriojYuvQz4AXCnmT0PfBN4DfAmYJZq2pabt5nZWqrVqT4p6cYW2kziAuRvBr7nisahjzrGcq5ImgLOAQ+4olng1WZ2LfBp4DuScq5JWXcMxqGP2+i9abetiyXDuA3z2FK3Jb2Iyig/YGY/BDCz02b2HzP7L/AtFmFq2A8ze879PQM85No87afp7u+Z3HI41gP7zOy0k6l1fVDf99bPFUkbgfcCHzbnVHWugz+57SepfLuvyyVDwzFoVR+qUok/AOwOZGtVF0uJcRvmsaRuO1/ZvcBRM/tqUB76LG8h84p3kl4q6SK/TfXA6TC9KZwbgR/llCOgZ0TUtj4cdX2fBj7qojOuA/7qXR45kLQO+Cxws5n9PSi/UtUavEhaDVwDPJNRjrpjMA1sULUI+yonxy9zyQG8A3jK3AppTrZWdbGkGPfTR6qn7U9T3W2nWmrzeqpp30GqV+7ud3LcDxxy5dPA8sxyrKZ6sn4AOOL7D1wOPAocd38va0EnL6F6JfglQVlWfVDdBGaBs1QjwNvr+k41db/bnSeHgDdnlmOGyofrz4973L4fdMfqALAPeF9mOWqPATDl9HEMWJ9LBlf+beAT0b7ZdLHUPyUlu1AoFDrGuF0ZhUKhUIgohrlQKBQ6RjHMhUKh0DGKYS4UCoWOUQxzoVAodIximAuFQqFjFMNcKBQKHeN/S7lN2vIrnEkAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x22bcef6db38>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from __future__ import division\n",
    "from osgeo import gdal\n",
    "import matplotlib.pyplot as plt\n",
    "from sklearn import cluster\n",
    "\n",
    "# read the data\n",
    "dataset = gdal.Open('datas/band1.tif')\n",
    "data1 = dataset.GetRasterBand(1).ReadAsArray()\n",
    "dataset = None\n",
    "\n",
    "dataset = gdal.Open('datas/band2.tif')\n",
    "data2 = dataset.GetRasterBand(1).ReadAsArray()\n",
    "dataset = None\n",
    "\n",
    "dataset = gdal.Open('datas/band3.tif')\n",
    "data3 = dataset.GetRasterBand(1).ReadAsArray()\n",
    "dataset = None\n",
    "\n",
    "\n",
    "dataset = gdal.Open('datas/band4.tif')\n",
    "data4 = dataset.GetRasterBand(1).ReadAsArray()\n",
    "dataset = None\n",
    "\n",
    "data1 = 1.0*data1\n",
    "data2 = 1.0*data2\n",
    "data3 = 1.0*data3\n",
    "data4 = 1.0*data4\n",
    "\n",
    "L = 1\n",
    "C1 = 1\n",
    "G = 2.5\n",
    "\n",
    "evi = G*(data3-data2)/(data3+C1*data2+L)\n",
    "\n",
    "data = evi.reshape(150*200,1)\n",
    "k_means = cluster.KMeans(4)\n",
    "k_means.fit(data)\n",
    "foo = k_means.labels_.reshape(150, 200)\n",
    "\n",
    "plt.contourf(foo,3)\n",
    "plt.colorbar()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.5.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
